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Intermediate hadronic meson loop contributions to J/ij), ip' — > 777c (7?7c) are studied apart from 
the dominant Ml transitions in an effective Lagrangian approach. Due to the property of the unique 
antisymmetric tensor coupling in V — > VP, the hadronic loop transitions provide explicit corrections 
to the Ml transition amplitudes derived from the naive "quenched" cc transitions via the coupling 
form factors. This mechanism interfering with the Ml transition amplitudes naturally accounts for 
the deviations from the Godfrey-Isgur model predictions in J/ip and tp' — » 777,3. It also predicts a 
small branching ratio of ip 1 — -> 777^, which can be examined by experimental measurements at BES 
and CLEO-c. 



I. INTRODUCTION 



Charmonium spectrum and decays of charmonium states are an ideal place for studying the strong 
interaction dynamics in the interplay of perturbative and non-perturbative QCD regime. In the past 
decades there have been significant progresses on the measurement of charmonium spectrum and their 
decays, which provide important constraints on phenomenological approaches. 

As the first charmonium state discovered in the history, J/ip has been one of the most widely studied 
states in both experiment and theory. As a relatively heavier system compared with light qq mesons, 
the application of a nonrclativistic potential model (NR model) including color Coulomb plus linear 
scalar potential and spin-spin, spin-orbit interactions, has provided a reasonably good prescription for the 
charmonium spectrum. This success is a direct indication of the validity of the naive "quenched" cc quark 
model scenario as a leading approximation in many circumstances. A relativised version was developed 
by Godfrey and Isgur 0] (GI model), where a flavor-dependent potential and QCD-motivated running 
coupling are employed. In comparison with the nonrelativistic model, the GI model offers a reasonably 
good description of the spectrum and matrix elements of most of the u, d, s, c and b quarkonia 0, 0|- 

On the other hand, there also arise apparent deviations in the spectrum observables which give warnings 
to a simple qq treatment and more complicated mechanisms may play a role. As pointed out in Ref. [J], 
the importance of mixing between quark model qq states and two meson continua may produce signifi- 
cant effects in the spectrum observables. By including the meson loops, the quark model is practically 
"unquenched" . This immediately raises questions about the range of validity of the naive "quenched" qq 
quark model scenario, and the manifestations of the intermediate meson loops in charmonium spectrum 
and their decays. These issues become an interesting topic in the study of charmonium spectrum with 
high-statistic charmonium events from experiment. An example is the newly identified state A(3872) 
and a possible assignment for it as a mixture of cc and DD* @, @] , or open charm effects 0] • 

In the recent years, the intermediate meson loop is investigated in a lot of meson decay channels [8l l9l. Hoi. 
[ill. Il2l. IT3L \\a\ flBL fl6j] as one of the important non-perturbative transition mechanisms, or known as final 
state interactiions (FSI). In particular in the energy region of charmonium masses, with more and more 
data from Belle, BaBar, CLEO-c and BES, it is widely studied that intermediate meson loop may account 
for apparent OZI-rule violations [HI, [H, [H, 0j EH (3 via quark-hadron duality arguement [l7l [l8l [l9j. 

In this work, we shall study the radiative decays of J/ip and ip' into 777c and 777^,. In the naive qq 
scenario, this type of decays is dominantly via magnetic dipole (Ml) transitions which flip the quark 
spin. For J/ip — > 777c and ip' — > 777^, where the initial and final state cc arc in the same multiplct, the 
spatial wavefunction overlap is unity at leading order, while 7/7 — > 7?7 C will vanish due to the orthogonality 
between states of different multiplets. In this sense, the former decays are "allowed" while the latter is 
"hindered" . However, the inclusion of relativistic corrections from the quark spin-dependent potential 
will induce a nonvanishing overlap between states of different multiplets such that the decay of ip' — > jrjc 
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is possible [13, SH H2, H • 

Theoretical studies of the heavy quarkonium Ml transitions with relativistic corrections are various in 
the literature [H, [HI, [H, H3, HI, HI, H3, HH . Relativistic quark model calculations show that a proper 
choice of the Lorentz structure of the quark-antiquark interaction in a meson is crucial for explaining the 
J ' — > 7?7 C data [27j • Systematic investigation of the Ml transitions in the framework of nonrelativistic 
effective field of QCD has been reported in Ref. [3fJ, where relativistic corrections of relative order v 2 are 
included. For J/tp — > 7?7 C the authors found Tj/^^ 1Vc = (1.5 ± 1.0) keV, which is in a good agreement 
with the data, but with quite large estimated uncertainties from higher-order relativistic corrections. 
Taking into account the transitions of ip' — * jr]c, the overall results for the Ml transitions still turn 
out to be puzzling [H, [Hj]. As studied by Barnes, Godfrey and Swanson Q in the NR model and the 
relativied GI model, although the models provide an overall consistent description of most of the existing 
charmonium states, theoretical results for the Ml transition have significant discrepancies compared with 
the experimental data Q. For example, in both NR and GI model, the predicted partial decay widths for 
J/ip — > 777c are as large as about two times of the experimental value, while for ip' — > ~/i] c , the theoretical 
predictions are about one order of magnitude larger than the data [3]. For ifj' — > 7?i c , although the 
predicted partial decay widths 0.17-0.21 keV are smaller than the experimental upper limit (< 0.67 keV), 
it is possible that the Ml transition is very different from the experimental measurement. 

Therefore, it is likely that there exist additional mechanisms beyond the cc transitions. This consid- 
eration thus prompts us to explore possible sources which can contribute to the charmonium radiative 
decay and cause deviations from the NR and GI model predictions, among which the intermediate meson 
loop transitions could be a natural mechanism. 

As follows, we first brief the calculations from the NR and GI models for the Ml transitions, and then 
introduce the formalisms for the intermediate meson loop contributions in Section II. The results and 
discussions will be presented in Section III. 



II. Ml TRANSITION IN NR AND GI MODEL 

The detailed study of the Ml transition was given by Barnes et al. in Ref. and here we quote their 
standard formula to incorporate the intermediate meson loop contributions which is to be introduced 
later. 

In Ref. [H, the partial decay width via Ml transition is evaluated by 

rM!(, 2S+1 I, - n' 2S ' +1 L' r+7 ) = ^^S^s-^m^E^ , (1) 

where n and n' are the main quantum number of the initial and final state charmonium meson; S (5"), L 
(£/) and J (J 7 ) are the initial (final) state spin, orbital angular momentum and total angular momentum. 
E 1 and Ef denote the final state photon and meson energy, respectively, while Mi is the initial cc meson 
mass. \ipi) and \ipf) arc the spatial wavefunctions of the initial and final state cc mesons, respectively. 

In the GI model, phase space factor Ef /Mi is not included though it is close to unity in many considered 
cases. In both GI and NR model, a recoil factor jo(kr/2) is included. We quote the results from Ref. [J] 
for future comparison. 

In order to incorporate the intermediate meson loop contributions, we derive the effective VjP cou- 
plings due to the Ml transition from Eq. fT]) by defining 

M fi (Ml) = ^e^P^Pte" , (2) 

where Pi and P 1 are four-vector momentum of the initial meson and final state photon, respectively, 
and £i and e 7 are the corresponding polarization vectors. From Ref. [J], we know that these extracted 
effective gv-yP couplings for J/ip^ip' — > ^rjc apparently overestimate the experimental data. Thus, the 
introduction of the intermediate meson loop contributions, which unquench the naive cc configurations, 
is supposed to cancel the Ml transition amplitudes via destructive interferences. 
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III. INTERMEDIATE MESON LOOP CONTRIBUTIONS 



The inclusion of the intermediate meson loops in meson decays somehow "unquenches" the naive quark 
model. A full consideration of such an effect requires systematic coupled channel calculations for e.g. 
the charmonium mass spectrum [3~H |. An interesting feature arising from the low- lying charmonia, such 
as r\ c , r]' c , J/ip, and ip' } is that their masses are lower than the open charmed meson decay channels. 
As a consequence, the lowest open charmed meson decay channels are expected to be dominant if they 
indeed account for contributions beyond the Ml transitions. This scenario turns to be consistent with the 
break-down of local quark-hadron duality, where the leading contributions to the sum over all intermediate 
virtual states are from those having less virtualitics. 

It should be pointed that intermediate states involving flavor changes turn out to be strongly sup- 
pressed. One reason is because of the large virtualities involved. The other is because of the OZI rule 
suppressions. Therefore, intermediate state contributions such as pir etc, are negligibly small. 

Following the above consideration, we thus investigate DD(D*), DD*(D*) and DD*(D) loops as the 
major contributions to J/tp — > jt) c , and tp' — > 777c 777^ as illustrated in Fig. [TJ We stress that although 
some of the vertices in the loop may violate gauge invariance, such as J /ipDD, the overall antisymmetric 
property is retained for the loops. The loop contributions hence only provide corrections to the VVP 
coupling strength for the external fields, but not change their antisymmetric tensor structure, no matter 

V is a massive vector meson or photon. Apart from the transitions in Fig. [TJ the contact transitions in 
Fig. [2] will also contribute to the decay amplitude. We show that the processes of Fig.[5]are gauge invariant 
by themselves. In brief, due to the property of the antisymmetric tensor coupling of VVP, where both 

V and P are external fields here, hadronic loop corrections are guaranteed to be gauge invariant in this 
effective Lagrangian approach. 

The detailed formulation is given in the following subsections. 



The transition amplitude for an initial vector charmonium (J/ip or ip') decay into 7?y c or 777^ via 
DD(D*) can be expressed as follows: 



where gi, gi, and 53 are the coupling constants at the meson interaction vertices (see Fig. [TJ). The four 
vectors, Pi, P 7 , and Pt are the momenta for the initial vector, final state 7 and pseudoscalar meson, 
respectively, while four- vector momenta, p\, pi, and p^ are for the intermediate mesons, respectively, 
and a\ = p\ — mf , «2 = p| — m 2i and 0,3 = p\ — m| are the denominators of the propagators of these 
intermediate mesons. 

As being studied in Ref. [l6j], this loop diverges logarithmically. Thus, a form factor to suppress the 
divergence and take into account the momentum-dependence of the vertex couplings is included: 



A. Intermediate DD(D*)+c.c. loop 




(3) 



where the vertex functions arc 




(4) 




(5) 



where n = 1,2 correspond to monopolc and dipolc form factors, respectively. An empirical argument 
applied here is that in the P-wave V — > VP decay the form factor favors a dipole form. We hence deduce 
the loop transition amplitudes with a dipole form factor. 
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Substitute the vertex couplings of Eq. (j4|) into Eq. {3]), the integral has an expression: 

" ~ ./ (2*)* J- (pl-mlM-mlM-ml) • 



With a dipole form factor, we have 



M fi = %-e a(ill3 ,P?&P?4 , (7) 



where 



Mi 



Mfl fX—x 



dx / dy- 



m 2 Jo Jo (4tt) 5 
y(A 2 - mf) 



A(mi,m 3 ,A) 
log — - — 

A(mi,m 3 ,m 2 ) 



(8) 



A(mi,m 3 ,A) 

where the function A is defined as 

A(a, b, c) = —(Mf - M 2 f ){\ -x-y)x + M 2 f x 2 + a 2 (l - x - y) - (M 2 f - b 2 )x + yc 2 . (9) 

In the intermediate meson exchange loop, coupling g 2 can be determined via the experimental infor- 
mation for D*° -» D°j(D*° -» L>° 7 ), i.e. 

2 127TMj,. 

_ — — J-.D*°-kD°7, (tUj 

where r^o^^o-y == (38.1 ± 2.9)% x T to t is given by experiment We neglect the contributions from the 
charged meson exchange loop since T £,,±^£ ) ± 1 = (1-6 ± 0.4)% x 96 keV is about two orders of magnitude 
smaller than r^o^^o^,. 

For coupling constant gi, especially gj/^uo, there are several methods suggested in the literature 
including quark model using heavy quark effective theory approach [35[, QCD sum rule [H, SU(4) 
symmetry and vector meson dominance (VMD) model [3q |. They typically give a value of order of one 
for 

9j/tjjDD- in this work, we adopt gj/^DD — 7-20 which is consistent with the value from Ref. [2a] . 
For the gD*Dr/ c coupling, we assume 

gD*°D°r lc = 5j/i/>Z3°*D° ■ (11) 



B. Intermediate DD*(D*)+c.c. loop 



As shown by Fig. [TJb), the transition amplitude from the intermediate DD*(D*) + c.c. loop can be 
expressed the same form as Eq. except that the vertex functions change to 



( rp = pa P V v 

11 — Mi a Pw r i t iP3 t 3 J 



T 2 

T 3 



iai I ,n a ' F P F V ' 

jn^ot'P'ii'v'Vl £ 2 r 7 t 7 

iS-P nan ,, n n a " P^" T) M " 
M f ta P V v P2 fc 2 P3 



-P' 



where /i,2,3 are the coupling constants. With a dipole form factor the integration gives 

9b 



M 



where 



fi 



/1/2/3 f 1 f , 

— it / dx / d y 

m 2 M f J J 



M, 



£aP[i 



pot P pH~i> 
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.4 



D 



'(4^)2 ( A2 A^ 



(12) 



(13) 



(14) 
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with 

A = \[{l-x- Z -){M*-M})+xM}] , 

B = -|(M? - M 2 f )[z(l - x)(M? - 3M|) - (M? - M})] , 

Ai = -xz(M? - Mf) + zm\ + ym\ + xm\ + (1 - x - y - z)A 2 . (15) 

In the above equation the intermediate meson masses 7774,2,3 are from the DD*{D*) loop, which are 
different from those in Eq. ([7J- /i,2,3 denotes the corresponding vertex coupling constants. 

In the DD*{D*) + c.c. loop, the coupling constant gj/^jj'jj is related to gjj^ub by the relation of the 
heavy quark mass limit (35| : 

/M d , (16) 

where Mjj corresponds to the mass ratio of Md/Md*. Similarly, we have g^u*D = 9ip' ddi 

For ip' -»■ 77y c and 777^, we assume that su/ D £, = gj/^ DD , 9d*°d° Vc = gj/^D*D, and «to»°£>X = 9t(,>d*d> 
which are consistent with the 3 Pq model [4[. In Table U the values of the coupling constants are listed. 



C. Intermediate DD*(D)+c.c. loop 



The transition amplitude from the intermediate DD*(D) + c.c. loop can contribute via charged in- 
termediate meson exchange. Treating the intermediate mesons as fundamental degrees of freedom, we 
eventually neglect the contributions from the non-zero magnetic moments of the D mesons. The charge- 
neutral loop is thus suppressed due to the vanishing D°D°'y electric coupling. Therefore, we only consider 
the charged meson loop contributions as shown by Fig. [He). The transition amplitude can also be ex- 
pressed in a form as Eq. ([3]) with the vertex functions 



Mi £ af3^Pi" £ 'i P3 £ 3 



i/a(Pi - 



P2) 
-P2) 



c 7 ■ 
£3 



(17) 



where /{ 2 3 are the coupling constants and and T(jy^) is the form factor. With a dipole form factor the 
integration gives 



M fi ee j^e^pvePpts" , 



(18) 



where 



/1/2/3 



dx 



dy 



(4tt) 2 



y(A 2 - m 2 2 ) 



A(mi,m 3 , A) 



log 



A(777l,7773, A) 

A (mi, 7773,7772) 



(19) 



It is interesting to note that the integral of the DD*(D) + c.c. loop has a similar form as that of 
DD(D*). However, we expect that contributions from this loop integral will be relatively suppressed 
since coupling f 2 is taken as the unit charge e = (47ra e ) 1 / 2 . 



D. Contact diagrams 



The contact diagrams of Fig. [2ja) and (b) (as an example in J/ip — > 7?y c ) arise from gauging the strong 
J/ifj(ip')D*D and r] c (ri' c )D* D interaction Lagrangians containing derivatives. The general form of the 



G 



transition amplitude of Fig. [2fa) and (b) can be expressed as follows: 

D*pol 

where ^(pl) is the form factor as before, and Ti(i = 1,2) are the vertex functions. For Fig. HJa), the 
expressions of Ti(i = 1,2) are: 

1 — Mi ^olPilv l fc 7 fc i P 2 t 2 r % t i fc 7 fc 27 ) (21) 

T 2 = 2ih 2 P f • e 2 , y 1 

where hi t % represent the J/ip(ip')D*D and rj c (r]' c )D*D coupling constants, respectively, and their values 
have been given in Tabled] 

Using the Feynman parameter scheme, the amplitude for Fig. [Ha) can be reduced to 



2i hl h 2 e f d* P2 fctfrfltf + P^^(-Pf + Pj ^)\{ml - A 2 ) 2 



M fi = U £a P^ 



Mi " ap »"J (2tt)4 (p 2 -m 2 )(p 2 -m 2 )(p 2 -A 2 ) 2 

2 ihl h 2 e r dH ~ *PtYf t + Pt^i-Pjf + w-*Pfn>-**r )](ml _ A2)2 



Mi opp j (2tt) 4 (l 2 -A 2 ) i 
'' ' ^-e^P^Pts^ f dx dy-L,(-^ - * ) (22) 



with 



and 



Mi 1 1 1 1 J J "(4tt) 2V 3A2 6A 2 

^-e a ^ v P^P^i , (23) 



A 2 = a: 2 M 2 - xM 2 f + xm\ + ym% + (1 - x - y)A 2 , (24) 



(4tt) 21 3A 2 6A 2 



9d = hih 2 e I dx / dyjj^ (—3 - ^-). (25) 



For Fig. IHb) , the vertex functions are 

T 

T 2 = 2e/i 2 e 7 • e 2 
The amplitude can then be reduced to 



71 = i^ir „ poc' 3 /? 1 ' 



(26) 



_ -2ihih2e b f d A p 2 p| 

M/l " Mi ^ £ -> J (2tt)4 (pi - m\)(pl - m 2 )(p 2 - A 2 ) 2 

-2ih x h 2 e DQ - „ /■ d 4 Z (i - sPi)" 



with 



Mi £ «/3^A%£ 7 y ( 27r )4( ? 2_ A3) 4' ( 28 ) 



A 3 = x 2 M 2 - zM 2 + im] + ym 2 + (1 - a; - y)A 2 . (29) 



Note that the integrand has an odd power of the internal momentum, the amplitude will vanish and has 
no contribution to the VVP coupling. 

The above deduction shows that only Fig. [3a) has nonvanishing contributions to the transition am- 
plitude. Meanwhile, gauge invariance is also guaranteed for the contact diagrams. The divergence of the 
loop integral is eliminated by adding the dipole form factor as in Fig. [1] 
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IV. RESULTS AND DISCUSSIONS 

Proceed to numerical results from the intermediate meson exchange loops, the undetermined quantities 
include the cut-off energy A in the dipole form factor and the relative phases among those amplitudes. 
The transition amplitude accommodating the Ml and intermediate meson exchange loops, i.e. DD(D*), 
DD*(D*), DD*(D), and the contact term, can then be expressed as 

Mfi = W [9VjP + 5 ° e<tfa + ~ 9he%5h + ^ + ade^e^PfsPpfe" , (30) 

where gv-yP is a real number and fixed to be positive. Couplings g a , g c and §d, calculated by the loop 
integrals can be complex numbers in principle. In this interested case, since the decay threshold of the 
intermediate mesons are above the initial meson (J/ip and if)') masses, the absorptive part of the loop 
integrals vanishes as a consequence. However, there might exist relative phases among those transition 
amplitudes. We hence include possible relative phases 5 a jj. c .d in the above expression. 

Note that the loop contributions are supposed to provide cancellations to the Ml amplitude which is 
real. We thus simply take S a ,b,c,d = or 7r. In this way, we have several phase combinations which are to 
be examined in the numerical calculation. Yet there is still a free parameter A to be constrained. 

We find that a reasonable constraint on the model can be achieved by requiring a satisfactory of the 
following conditions: i) For either constructive (5 — 0) or destructive phases (5 — tt), the same value 
of A is needed to account for J/ip — > "frj c , ip' — > 7?7 C simultaneously, ii) The value of A is within the 
commonly accepted region, 1.5 ~ 2.5 GeV. hi) The prediction for ip' — > 777^, with the same A is well below 
the experimental upper limit, BR(ip' — > 777^) < 2.0 x 10~ 3 [3]. 

Imposing the above conditions on fitting the A parameter, we obtain A = 2.39 GeV as the best fit 
with 5 a ,b,c,d = t, i-c contributions from the loop integrals provide cancellations to the Ml transition 
amplitudes and there is no need for abnormal relative phases among the intermediate meson exchanges. 

The numerical results for the intermediate meson exchanges have some predominant features. We find 
that the DD(D*) and DD*(D*) loops have relatively large contributions while the DD*(D) loop is quite 
small. The contributions from the contact term are negligibly small in J/ip — > jr) c and ip' — > 777^,, while 
relatively large in ip' — > jTjc- 

In Table Hi] the fitted branching ratios are listed and compared with the GI model Ml transitions. We 
also list the exclusive contributions from the triangle diagrams of Fig. [1] and contact diagrams of Fig. [5] 
as a comparison. For J/ip — > 777c we find that the magnitude of the meson loop amplitude is smaller 
than the Ml amplitude, while for ip' — + r )-q Cl the absolute loop amplitude turns to be larger than the Ml. 
With A = 2.39 GeV, we obtain T(J/ip — » 7?y c ) = 1.59 keV which is located at the upper limit of the 
experimental data, T exp (J/ip — > 7?/ c ) = (1.21 ±0.37) keV Q. For ip' — > 777c we have T(il>' — > 7?/ c ) = 0.86 
keV, which is agree well with the data, T exp (ip' -> 777J = (0.88 ± 0.13) keV 0]. Taking into account 
the still-large uncertainties with the data for J/ip — > 7?y c , the inclusion of the intermediate meson loop 
contributions significantly improves the theoretical results. 

With the fixed A, the partial decay width for ip' — > 777^. is calculated as a prediction. The pure Ml 
transition predicts r^/ 1 (-0' — > 77^,) ~ 0.17 keV while the hadronic loops contribute Thl(^' — * IVc) — 
0.054 keV. The cancellation from the hadron loops thus leads to r a u(ip' — > 77^,) ~ 0.032 keV which is well 
below the experimental upper limit, 0.67 keV 0|. Note that in all these three channels, the hadronic loop 
cancellations from the real part of the amplitudes possess the same relative sign to the Ml amplitudes. 
This makes the decay of ip' — > 777^, extremely interesting. As the pure Ml transition still predicts a 
sizeable partial width about 0.17 keV while our hadronic loop cancellation predicts a much smaller value, 
improved measurement of this quantity will help us gain further insights into the decay mechanisms. 

For other relative phases, we find that there does not exist a common value for A to fit the data for 
J I ip — * 7?7 C and ip' — > 7?y c simultaneously. 

To summarize, in this work we have studied the hadronic meson loop contributions to the J/ip and ip' 
radiative decays into 77y c or 777^.. In the framework of effective Lagrangian phenomenology, the interme- 
diate meson exchange loops provide corrections to the leading couplings extracted from potential quark 
models. In comparison with the NR and GI model, the meson loop contributions turn to cancel the NR 
and GI amplitudes. It is interesting to see that the meson loop contributions in J/ip — > "/rj c is smaller 
than the Ml transition in magnitude, while in ip' — > 7?y c the situation is opposite. Note that the pure Ml 



s 



contribution in ip — > 777c is about one order of magnitude larger than the experimental data, the meson 
loop contributions turn out to be even larger. This mechanism suggests significant cancellations between 
the Ml and meson loop amplitudes. It raises questions on the naive qq solution for the meson spectrum, 
and could be a manifestation of the limit of the quenched quark model scenario. 

As a prediction from this model, we calculate the partial decay width of ip' — > ji]' c . It gives a value 
about one order of magnitude smaller than the experimental upper limit. Improved measurement of this 
decay channel is strongly recommended. 

It is interesting to note that our model results are similar to those from a relativistic quark model 
calculation by Ebert, Faustov and Galkin [27j . who find that a proper choice of the Lorentz structure 
of the quark-antiquark interaction in a meson is crucial for accounting for the Ml transition data. In 
our approach we extract the effective couplings from the NR and relativised GI model and then combine 
it with the gauge invariant meson loop corrections. The validity of this approach is guaranteed by the 
property of the unique antisymmetric tensor coupling for VVP fields. In the framework of effective 
Lagrangian phenomenology the corrections to the leading contributions are introduced as coupling form 
factors. 

Although we also observe strong sensitivities of the hadronic loop contributions to the cut-off energy 
A, the advantage of this approach is that the number of parameters is limited. In fact, there is little 
freedom for the effective couplings at vertices. By a coherent study of J /ip and tp' — » 7?y c , we find that the 
constraint on the A value is very tight. Certainly, it should be noted that our treatment of the relative 
phases is empirical though the favor of a destructive phase between the Ml transition amplitudes and 
hadron loops turns to be consistent with what one naturally expects. Note that it has been shown in 
Ref. [23[ that a proper modification of the color Coulomb potential strength will simultaneously account 
for the branching ratios for J/ip — > 777c and tp' — > 7?7 C . It seems to support that the intermediate hadronic 
meson loops are responsible, at least partly, for such a modification, and hence break down the naive qq 
scenario. 

The study of non-perturbative effects arising from intermediate meson loops in heavy quarkonium de- 
cays has attracted a lot of attention recently. Although such approaches still experience large uncertainties 
from the divergent behavior of the loop integrals, we expect that improved experimental measurements 
with high statistics, such as at BES and CLEO-c, provide more and more stringent constraints on the 
hadronic loops. Thus, insights into the effective degrees of freedom within hadrons and their decay mech- 
anisms can be gained. This requires systematic analysis of both spectroscopy and coupled channels for 
which more and more theoretical efforts are undergoing. 
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Coupling constants 


\9.J/i>DD 




\9j/tPD*d\ 


\9d*Dt) c 1 


\9tP'd*d\ 


\9D*Dr)' a 1 


\gD'D~,\ 


numerical value 


7.20 


7.20 


4.34 


4.34 


3.64 


3.64 


6.86 



TABLE I: The absolute values of coupling constants for the effective vertex interactions. Their relative phases 
are determined by the SU(4) flavor symmetry. 




(a) 




(b) 




FIG. 1: Schematic diagrams for J/ip — ► ^r\ c via (a) DD(D*), (b) DD*(D*) and (c) DD*(D) intermediate meson 
loops. Similar processes occur in ip' — > ^r\ c and ^r\ a . 
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Initial meson 


J/V(l 3 5i) 


V'(2 3 5i) 


Final meson 


ijc(l l So) 


??c(2 1 So) 




r£? (keV) 


2.9 


0.21 


9.7 


T%\ (keV) 


2.4 


0.17 


9.6 


T C (keV) 


~ 


~ 


0.04 


IYh (keV) 


0.096 


0.063 


17.91 


Thl (keV) 


0.083 


0.054 


16.20 


T aii (keV) 


1.59 


0.032 


0.86 


T exp (keV) 


1.21 ±0.37 


< 0.67 


0.88 ±0.13 



TABLE II: Radiative partial decay widths given by different processes are listed: F M j and T < £ 11 are the Ml 
transitions in the NR and GI model, respectively Q; Tth are inclusive contributions from the triangle diagrams 
(Fig-HJ; Tc are from contact diagrams (Fig. [2}; Thl denote the inclusive contributions from all the intermediate 
hadronic loops; while F a ;; are coherent results including the Ml in the GI model and intermediate hadronic loops. 
The experimental data are from PDG2006 [3J. The results are obtained at the cut-off energy A = 2.39 GeV. 





(b) 



FIG. 2: The contact diagrams considered in J/ip — > 7^. Similar diagrams are also considered in ip' — » 'yr] c (ri' c ) 



